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ABSTRACT 



We describe radio continuum and spectral line observations of the source G320.4- 
01.2 (MSH 15-52) and the coincident young pulsar B1509-58 made with the Australia 
Telescope Compact Array. Based on a comparison between X-ray and radio observa- 
' tions, we argue that the two main radio components of G320.4-01.2 are a single super- 

\Q | nova remnant (SNR), which Hi absorption indicates is at a distance of 5.2 ± 1.4 kpc. 

£Nl ■ A high-resolution correspondence between radio and X-rays argues that the pulsar 

is interacting with the SNR via an opposed pair of collimated outflows. The outflow 
itself is seen as an elongated X-ray feature surrounded by a highly polarized radio 
On . sheath, while the interaction with the SNR manifests itself as a ring of radio/X-ray 

knots within the optical nebula RCW 89. We reject the hypothesis that the pulsar 
. outflow powers the entire RCW 89 region. 

SNR G320.4-01.2 and PSR B1509-58 agree in distance and in rotation measure, 
and appear to be interacting. We thus conclude that the two objects are associated 
and have an age of ^1700 yr. We propose that the SNR resulted from a high-energy or 
low-mass supernova which occurred near the edge of an elongated cavity. Such a model 
\ can account for the SNR's bilateral appearance, its large apparent age, the significant 

offset of the pulsar from the SNR's centre and the faintness of the pulsar-powered 
nebula at radio wavelengths. 

Key words: ISM: individual: G320.4-01.2, RCW 89 - ISM: jets and outflows - 
■ pulsars: individual: B1509-58 - radio continuum: ISM - shock waves 



1 INTRODUCTION 

A massive star ends its life in a supernova. This produces 
a supernova remnant (SNR) and often also an associated 
neutron star, the latter sometimes observable as a pulsar. 
However, associations between SNRs and pulsars are rare: 30 
years of effort have provided fewer than ten convincing cases 
(e.g. Kaspi 1996). Such associations can clarify questions 
regarding pulsar velocities, ages, magnetic fields and initial 
spin periods, and can help us to understand the evolution 
and appearance of SNRs. Thus establishing the validity of 
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an association is of great interest, as are subsequent studies 
of a particular pulsar/remnant pairing. 

Of such systems, one of the best studied but least un- 
derstood is that involving the peculiar SNR G320.4-01.2 
and the energetic young pulsar B1509-58. In this paper we 
present a new radio study of this intriguing pair of objects. 



1.1 Previous Observations 

1.1.1 Radio 

The radio source G320.4-01.2 (MSH 15-52, Kes 23) w as de- 



tected in the ear l iest radio surveys of the southern sky ( Mills 
Sice fe Hill 1960|: IfComesaroff 196&! ; |fflll 1968] ; |Kesteven 1961 
Day, Thomas fc Goss 1969| ), and was soon classified as a 



SNR on the basis of its non-thermal spectrum and lack 
of hydrogen recombination lines (Milne et al. 1969; Milne 
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Figure 1. A radio/X-ray comparison of 0320.4-01 2. The grey- 



scale c orres ponds to 36-cm MOST observations (Whiteoak & 
Green l [j96[ ), while the contours represent smoothed ROSAT Po- 
sition Sensitive Proportional Counter (PSPC) data of Trussoni 



et al. (L996) [their Figure 2(a)]. Contour levels (in arbitrary units) 
are at levels ot O.b. 1. l.b. 'I. b. 10. 2 0. 30 and 40. Selected regions 
designated by Trussoni et al. (1996) are indicated. Here (and in 
subsequent figures), the position of PSR B1509-58 is marked by 
a "+" symbol. G320.6-01.6 is to the south-east of G320.4-01.2 
and is outside the field of view shown here. 



1970; shaver fc Goss 1970c ) . Early images of the source 



flShaver fc Goss 1970b| ; [Milne 1972 



Milne & Dickel 197E 



showed (at least) two distinct components: a bright double 
source to the north-west, and fainter extended structure to 
the sout h-east. Subsequent observations at s ub-arcmin res- 

[Manchester & 

Durdin 



olution (|Caswcll, Milne fc Wellington 198l[ 

1983| ; |Whiteoak fc Green 1996] ) showed the north- 



western component to be dominated by a central clumpy 
ring, and the south-eastern source to be comprised of com- 
plex filamentary structure, as shown in Figure |l|. Both com- 
ponents are extended along a direction parallel to the Galac- 
tic Plane, with a gap between them. 

The radio morphology of G320.4-01.2 is unusual. Its 
two well-separated components can only be considered a 
single shell in the broadest of interpretations; the unusual 
collections of knots and filaments, and the distinct gaps in 
emission to the north-east and south-west, need to be ex- 
plained. The regio n is probably the site of many super novae 
in the recent past (Lortet, Georgelin & Georgelin 1987), and 
one possibility is that the components of G320.4-01.2 may 
actuall y be multiple SNRs. Wide-field observations of the 
regi on ( Milne, Caswell & Haynes 1993; Whiteoak & Green 



1996) rsveal a pair of polarized radio arcs ~30 arcmin to 
the south-east of G320.4-01.2, which have been designated 
G320.6-01.6. It is not clear whether this structure is another 
part of G320.4-01.2 or is a separate SNR. 




2 X-ray 

The first X-ray image of the region revea led a surprise: a 



150-ms pulsa r B1509-58 within the SNR (Beward fc Harn- 



den Jr. 1982 ). The pulsar was subsequently detected in the 



radio ( Manchester, Tuohy fc D'Amico 1982 ) and 7-ray ( Wil- 
son et al. 1992;; |Ulmer et al. 1993|) bands. Radio timing ob- 



servations ( [Manchester, Durdin fc Newton 1985 ; Kaspi et al 



1994) show PSR B1509-58 to have the highest period deriva- 
tive (P = 1.5 x 10~ 12 ), the second lowest characteristic age 
(t c = 1700 yr) and the third highest spin-down luminosity 
(E = 1.8 x 10 37 erg s -1 ) of any known pulsar. 

The X-ray morphology of the SNR is shown in Fig- 
ure 0. The pulsar is embedded in a 5-10 arcmin non- 
thermal nebula (region "SN" o f Figure fl|) (Reward et al 



198J; preiveldinger et al. 1995| ; [Marsden et al. 1997] ), pre- 



sumed to be synchrotron emission from a pulsar wind nebula 
(PWN). There is no obvious radio counterpart to this PWN. 
To the north-west of the pulsar is an extended thermal X- 
ray s o urce (region "NN") (peward et al. 1983 ; Trussoni et al 



Tamura et al. 1996), whose position and m orphology 



closely corre spond to the brightest radio emission ( Brazier & 



Becker 1997). A diffuse bridge of emission joins the "SN" and 



"NN" regions. Another bridge connects the "SN" region to 
X-ray emission south-east of the pulsar (region "T" ) . Region 
"T" appears to coincide in position and morp hology with the 



south -eastern component of the radio SNR ( Trussoni et al 



1996). The lowest level X-ray contour envelopes the entire 
SNR (region "CDN"). 



1.1.3 Optical and Infrared 

The bright north-western X-ray and radio source coincides 
with the collection of irregular filaments which make up the 
Ha nebula RCW 8 9 ( |Rodgers, Campbell fc Whiteoak 196c| ; 
Seward et al. 1985 ). O ther nearby nebulosities are probabl y 
unrelated to the SNR ( Lortet, Georgelin fc Georgelin 1987 ) 



The infrared point source IRAS 15099-585 6 is coincide nt 
with the pulsar, but its nature is uncertain (Arendt 1991). 



1.2 Are G320.4-01.2 and B1509-58 associated? 

PSR B1509-58 is a young pulsar located near the cen- 
tre of a SNR. At first glance, a physical association with 
G320.4-01.2 seems assured. However, such a claim should be 
treated with caution: a young pulsar need not have an asso - 
ciated SNR (|Kafatos et al. 1980| [Braun, Goss fc Lyne 1989|), 



and line-of-sight coincidences are rife (Gaensler & Johnston 
1995a,b). There are several criteria by which an association 
can be judged. The most useful are agreement in distance 
and age estimates, a reasonable implied transverse velocity 
for the pulsar, evidence from the pulsar's proper motion, 
and any indication of an interaction between the pulsar and 
the SNR. 



Hi absorption towards RCW 89 (Caswell et al. 197E 



puts the SNR at a distance of 4.2 kpc, slightly less than 
the pulsar's dista nce of 5.9=1=0.6 kpc as de rived from its dis- 
persion measure ( Taylor fc Cordes 1993 ). A more serious 



discrepancy lies in the ages of the two systems: unless one 
invokes a time- variable br aking index, PSR B 1509-58 can 
be no older than 1700 yr (Kaspi et al. 1994). But assuming 



© 1998 RAS, MNRAS 000, 



Radio observations of SNR G3204-01.2 and PSR B1509-58 3 



standard supernova and i nterstellar medium (ISM) param- 
eters, Seward et al. (1983) calculate an age for the SNR in 



the range (6 — 21) x 10 yr. Proper motion measurements 
of the optical filaments in RCW 89 also imply a remnant 
age mu ch grea ter than that of the pulsar ( van den Bergh & 
Kampe' 1984). One way of reconciling the SNR and pulsar 



ages is by invoking rapid expansion of the SNR shell into 
a cavity, followed by an encounter with dense surrounding 
material (Seward et al. 1983; Manchester 1992). 



Calculating the transverse pulsar velocity implied by an 
association requires a determination of the pulsar birthplace. 
This is straightforward when the SNR is a near-circular shell, 
but is more difficult for the complex morphology of G320.4- 
01.2. Suffice to say that the pulsar's current position well 
within the SNR's perimeter implies a projected velocity less 
than~3000 kms~\ No proper motion has yet been detected 



flKaspi et al. 1994[ ). 

In order to expl ain R CW 89's bright appearance in X- 
rays, Seward et al. (1983) suggested that it is power ed by 
an outflow of energy from the pulsar. Tamura et al. (199C) 
revive this interpretation, showing that the X-ray bridge 
joining the pulsar and RCW 89 has a non-thermal spec- 
trum. They propose that this bridge is a jet flowing out from 
t he pu lsar and illuminati ng RC W 89. Manchester & Durdin 
(1983) and Manchester (1987) describe a related idea, not- 



ing that the two components of the radio remnant can be 
mapped onto opposed cones emanating from PSR B1509- 
58. They suggest that the radio morphology of the SNR may 
be generated by opposed beams of particles originating from 
the pulsar and colliding with t he surface of a surrounding 
cavity. Brazier & Becker (1997) provide X-ray evidence for 



this model, pointing out that the PWN has a cross-shaped 
appearance, possibly corresponding to a highly inclined ver- 
sion of the "torus + jets" geometry of the Crab Nebula in 
X-rays (Hester et al. 1995). The "jets" align with the axis of 



proposed outflow (presumed also to be the pulsar's rotation 
axis) and one of them maps onto the ring of X-ray knots 
seen at the peak of RCW 89. 

There are two main problems with the argument that 
the RCW 89 region is powered purely by an outflow from 
the pulsar. First, the total thermal energy in RCW 89 
(4 x 10 49 erg; Tamura et al. 1996) exceeds the entire en- 
ergy lost by the pulsar over its lifetime for any but very 
rapid initial periods. Second, the mechanism by which a 
large thermal nebula might be created at the te rmination of 
such an outflow is unclear (Tamura et al. 1996). Both these 
issues must be resolved if this model is to be considered a 
possibility. 



2 NEW OBSERVATIONS 

Clearly, there are many issues still to be resolved involving 
SNR G320.4-01.2 and PSR B1509-58. This has prompted us 
to carry out a set of new radio observations, using the Aus- 
tralia Telescope Compact Array (ATCA; Frater, Brooks & 
Whiteoak 1992) , a synthesis telescope located near Narrabri, 
NSW, Australia. The ATCA consists of five moveable 22-m 
antennas on a 3-km east-west track, with a fixed sixth an- 
tenna stationed 3 km further west. Five types of observations 
were carried out, as detailed in Table |l|: 



(i) continuum observations at 20 cm (1.3 GHz) and 6 cm 
(4.8 and 5.8 GHz); 

(ii) pulsar-gated observations; 

(iii) observations in the H I spectral line; 

(iv) observations in the 1720-MHz OH maser line; 

(v) observations in the H140a recombination line 
(2371.1 MHz). 

H I and ungated continuum observations at 20 cm con- 
sisted of two fields covering the entire SNR. Continuum ob- 
servations at 6 cm were a mosaic of eight fields, covering 
RCW 89, a region around B1509-58 and part of the south- 
eastern component. OH, H140« and pulsar-gated observa- 
tions consisted of a single pointing towards RCW 89. All 
four Stokes parameters were recorded in continuum, while 
total-intensity alone was recorded for line observations. 
The flux density scale of all data was tied to the revised 
scale of Reynolds (1994) using observations of PKS B1934- 
638. Antenna gains were determined using regular observa- 
tions of MRC 1613-586 (20-cm continuum, Hi and H140a), 
PKS B1740-517 (6-cm continuum) and PKS B0823-500 
(OH). 



3 DATA REDUCTION 

Reductio n and analysis were carried out within the MIRIAD 
package (3ault, Teuben & Wright 1995). Data were edited 



and calibrated using standard techniques (Sault & Killeen 



199S 



The field of the calibrator MRC 1613-586 contains 
several weak confusing sources, and antenna gains for this 
source were determined using a model consisting of both the 
calibrator and these weaker sources (cf. Gaensler, Manch- 
ester & Green 1998a). 



3.1 Continuum data 

At both 20 and 6 cm, a mosaic im age of the field was formed 
using multi-frequency synthesis (Sault & Wieringa 1994) 



and uniform weighting. Data from the 6-km antenna were 
excluded. Each ima ge was deconvolved us ing the maximum 
entropy algorithm (pull fc Daniell 1978), all pointings in 



the mosaic being han dled simultaenously (sault, Staveley- 
Smith & Brouw 1996). The resulting models were smoothed 
using a Gaussian restoring beam, and then corrected for the 
mean primary response of the ATCA antennas and for the 
appropriate mosaic pattern. The resolution and noise in the 
final images are given in Table |^. 

The pulsar-gated data were analyse d by de-dispersin g 



at a dispersion measure of 253.2 pc cm (Kaspi et al. 1994) 



On- and off-pulse datasets were formed, and the off-pulse 
data subtracted from the on-pulse data so that only pulsed 
emission remained. 

Images of the region were also formed in Stokes Q, U 
and V. At 20 cm, Faraday rotation across the observing band 
is significant (~4 radians across 128 MHz in some regions), 
and can cause depolarization of the emission. This effect 
was minimised by making multiple pairs of Q and U images 
across the 20-cm band, each of bandwidth 8 MHz. At 6-cm 
bandwidth depolarization is minimal (< 0.5 radians), and 
we made a single pair of Q and U images for the entire 
frequency band. 
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Table 1. ATCA observations of G320.4-01.2. Except where noted, all observations were in continuum mode, in- 
volving 32 channels across a 128 MHz bandwidth. 



Date 


Array 


Maximum 




^2 


Time on 




Config 


Baseline (m) 


(MHz) 


(MHz) 


Source (h) 


iyyo inov uo 




oyoy 


1 'XAA 




i 


1995 Nov 11 


6A 




1721 6 




3 


1996 Feb 01 


0.75B 


765 


1344 


1420 a 


3 


1996 May 10 


1.5D 


1439 


1376 


1420 a 


16 


1996 May 12 


1.5D 




4790 


5824 


6 


1996 May 13 


0.75D 


719 


4790 


5824 


9 


1996 May 16 


0.75D 




1344 


1420 a 


12 


1996 May 17 


0.75D 




4790 


5824 


12 


1996 Nov 10 


0.75A 


735 


2371 c 


1376 d 


13 


1997 Apr 18 


0.375 


459 


1344 


1420 a 


11 


1997 Apr 20 


0.375 




4800 


5824 


13 



a Hi observations (1024 channels across 4 MHz) 

b OH observations (1024 channels across 4 MHz); these observations also formed part of the survey of Green et al. 

(1997), but have been reprocessed here 

c H140c? observations (256 channels across 8 MHz) 

d Pulsar-gated observations (similar to continuum mode, but where data were sampled every 4.7 ms and then folded 
at the apparent period of PSR B1509-58) 



Although mosaic images can be formed in polarization 
as described for total-intensity data above, Q, U and V data 
are not positive everywhere and thus cannot be deconvolved 
using maximum entropy techniques as used for Stokes I. 
For our 20-cm data we thus deconvolved images of polar- 
ized emission using the "individual" approach (see Sault & 
Killeen 1998), whereby each pointin g is deconvo lved sepa- 
rately using the CLEAN algorithm (Clark 1980), and final 
images in Q, U or V are then formed by combining each 
CLEANed pointing. However the "individual" approach re- 
sults in an image greatly lacking in extended str ucture 



when compared to "joint" deconvolution techniques (Corn- 
well 19|i|). At 20 cm most of the emission is well sampled 
by the unmosaiced u — v coverage, and this effect is mini- 
mal. However at 6 cm the image quality resulting from the 
individual approach is poor. Thus at 6 cm we formed dirty 
polarization images using the joint algorithm as in total in- 
tensity, but made no attempt to deconvolve them. 

A linear polarization image L was formed from each pair 
of Q and U images and corrected for non-Gaussian noise 
statistics. At 20 cm, each L image was clipped where the 
polarized emission or the total intensity was less than 5a, 
and the final image formed by taking a mean across the 
observing band. 



3.2 Spectral Index 

Existing spectral index studies of G320.4-01.2 have had low 
spatial resolution and have produced a wide range of con- 
flicting results (see du Plessis et al. 1995). But using our 
high-resolution 20- and 6-cm observations, we can accurately 
calculate the spatial distribution of a (where S„ cx v a ) for 
G320.4-01.2. 

When computing spectral indices, it is crucial that the 
images being compared contain the same spatial scales. The 
u — v coverage of our two datasets is quite different, and 
spatial filtering must be applied to each so that their spatial 
scales match. Assuming that the deconvolution process at 



each frequency accurately fills in the u — v plane in some 
range w mllx < (u 2 + u 2 ) 1//2 < u> max , this filtering can be sim- 
ply achieved by smoothing the higher resolution image to 
the resolution of the lower, and discarding extended struc- 
ture in whichever image extends to shorter projected base- 
lines. The spatial distribution of spectral index can then 
be determined usin g the method of "T-T" (temperature- 
temperature) plots (Turtle et al. 1962). In this method, two 
images are plotted against each other on a pixel-by-pixel 
basis in some small sub-region of the image. Provided any 
DC offset is constant over this region, the data will follow a 
linear relationship , whose slope m is related to the spectral 
index, a, by 



lnra 



ln(^i/V 2 ) ' 



(1) 



A linear fit to the data gives a value of m and thus of a. 

We applied this method to our 20- and 6-cm data on 
G320.4-01.2 to produce a pair of images of resolution 24 arc- 
sec. Because there are many pixels per beam area in each 
image, adjacent pixels are correlated, which will cause the 
uncerta inty in the l inear fit to a T-T plot to be underesti- 
mated (Green 199C). To correct for this effect, each image 
was re-binned by choosing every ith pixel, where there are 
i pixels per beam in each dimension. T-T plots of v\ vs V2 
were then produced in square regions of side 2 arcmin, and 
a linear least squares fit applied to each plot to determine 
m, and hence a. 



3.3 Spectral line data 

For the H I observations, continuum emission was subtracted 



from the line data in th e u — v plane (van Langevelde & 



Cotton 1990; Sault 1994), and spectra were then smoothed 
to a velocity resolution of 3.3 km s . Diffuse emission was 
removed by discarding u — v spacings shorter than 1 kA, 
and a line cube was then formed for velocities between -200 
and +200 km s" 1 (LSR). The peak emission in the resulting 
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images was deemed sufficiently faint that deconvolution was 
not required. 

The H I cube was then weighted by multiplying it by the 
20-cm continuum image. Hi spectra towards sources of in- 
terest were generated by integrating over the relevant region 
of the cube, and then renormalised to give units of fractional 
absorption. 

The rms noise, a, in each spectrum was determined from 
the flux in line-free channels. The presence of extended H I 
emission in the Galactic Plane raises the thermal noise in 
the Hi line to approximately double this value (e.g. Dickey 
1997), and so we choose 6cr as the threshold above which 
features are considered believable. 

Continuum emission was subtracted from OH and 
H140q data as for Hi. Cubes were formed at various ve- 
locity resolutions and then searched for emission. 



4 RESULTS 

4.1 Continuum Images 

Total-intensity images of SNR G320.4-01.2 at 20 and 6 cm 
are shown in Figures |^ and ^ respectively. The brightest re- 
gion is that coincident with RCW 89, and can be seen to 
form a ring of knots (region 4 in Figure ^) superimposed on 
a horseshoe-shaped plateau of emission (region 5). A bridge 
of emission running parallel to the Galactic Plane connects 
RCW 89 with another clump of emission, comprised of a 
compact source at its south-west end (region 3) and becom- 
ing more diffuse to the north-east (region 2). South-east of 
RCW 89 is a twisted, filamentary arc (incorporating regions 
7 and 9), with one distinctly compact region near its north- 
ern end (region 8). 

The 20-cm flux density of G320.4-01.2 and its associ- 
ated uncertainty are given in Table |^, and were measured by 
integrating polygons around each component of the source 
then applying a background correction. The shortest an- 
tenna spacing at 20 cm is 31 m, corresponding to a spa- 
tial scale of ~25 arcmin. However, the mosaic allows us to 
recover information on larger scales (Cornwell 1988; Sault 
et al. 1996), in this case up to ~35 arcmin. If the SNR 
contains structure larger than this, these missing spacings 
will result in a significant underestimate in the true flux 
density. However, the 20-cm flux densities in Table |^ are in 
excellent agreement with s i ngle-dish ob servations at this fre- 
quency flMilne et al. 1969| ; |Milne 1972j ). Furthermore, there 
is no negative bowl observable around the SNR, as would 
be expected if large-scale structure were missing. Finally, 
our ATCA image, when smoothed to the appropriate reso- 
lution, is indistinguishable from a 20-cm image made with 
the Parkes 64-m radio telescope (A. R. Duncan 1996, private 
communication). All these points are good evidence that the 
20 cm ATCA observations contain all the flux and structure 
from the SNR. 

At 6 cm, we are sensitive only to scales smaller than 
9 arcmin, even after mosaicing. Thus we expect the flux 
density determined in Table ^| (for the northern component 
only) to be an underestimate. Indeed if we take the flux 
density at 20 cm and scale it by a mean spectral index a — 
—0.45 (Milne et al. 1993), we expect a 6-cm flux density 
for the northern region of 15 Jy, substantially greater than 



Table 2. Observational and derived parameters for ATCA con- 
tinuum data. 



Wavelength (cm) 



20 



Resolution 

Largest angular scale 

in image 
Measured rms noise 

(|ijy beam ) a 

G320.4-01.2 : 

Flux density (Jy) a 



PSR B 1509-58 : 

Flux density (mjy)° 
Position (J2000) 



24'.' 1 X 2(y.'8 
35' 

250 
75 



28 ±1 
18 ±1 
46 ±2 



(Stokes I) 
(Stokes V) 



(north) 
(south) 
(total) 



1.3 ±0.3 
RA 15 h 13 m 55?61 ±0?02, 
Dec -59°08'08'.'67±0"26 



15'.'0 x 10V4 
9' 

250 
100 



10 ± 1 



< 0.5 



G320.6-00.9 (Source 1 in Figure || : 

Flux density (Jy) a 0.35 ± 0.05 0.1 
Position (J2000) RA 15 h 15 m s , Dec -58°46'53" 

a 1 jansky (Jy) = 10" 2B W m" 2 Hz" 1 

the value in Table |j| This missing flux is clearly evident 
in Figure ^, which, compared to Figure ^, lacks significant 
extended structure. 

The position and 20-cm flux density of PSR B1509-58 
are given in Table |^ and were determined by fitting in the 
u — v plane to the pulsed dataset formed from the pulsar- 
gated data. This position is marked by a "+" symbol on 
appropriate Figures. The puls e-prof ile is similar to that ob- 
tained by Manchester et al. ( 1982 ), consisting of a single 
main pulse occupying ~15 per cent (FWHM) of the period. 
The pulsar was not detected in (ungated) 6-cm observations. 

The position of the pulsar given in Table jsj is consistent, 
within the uncertainties, with that derived by Manchester 
et al. (1985) 14.4 yr earlier. This allows us to put 1-cr upper 
limits on the pulsar's proper motion of 39 mas yr" 1 in RA 
and 52 mas yr" 1 in Dec. The corresponding 1-<t error ellipse 
for the distance travelled in 1700 yr is shown in Figure ^, 
limiting its birth-place to somewhere south of RCW 89 and 
well within the SNR. 

An image of the region surrounding the pulsar is given 
in Figure ^, where the grey-scale range has been chosen to 
show faint radio emission in the region. The pulsar is seen 
to lie within a tongue-like plateau of emission extending 6 
arcmin south-east from RCW 89 (region 6 of Figure ^]); a 
channel of depressed emission in this tongue can be seen sur- 
rounding and to the south-east of the pulsar. This tongue 
and the channel within it are clearly evident in MOST ob- 
servations of the same r egion (Manchester & Durdin 1982 ; 
Whiteoak & Green 1996), but under conditions of lower res- 
olution, poorer sensitivity and reduced dynamic range. It 
is also just visible in our 6-cm data, but is poorly imaged 
as a result of the lack of short spacings. X-ray emission in 
the region (as imaged by ROSAT PSPC) is shown as con- 
tours in Figure ^. A comparison between X-ray and radio 
data clearly shows that a narrow X-ray feature (part of the 
PWN powered by PSR B1509-58) runs along the channel of 
reduced radio emission. 

In Figure H we compare radio emission from the centre 
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15 h 16 m 15 m 14 m 13 m 12 m 15 h 16 m 15 m 14 m 13 m 12 1 

RA (J2000) RA (J2000) 



Figure 2. Grey-scale and contour images of SNR G320.4-01.2 at 20 cm, corrected for the primary beam response of 
the telescope. The grey-scale range is to 75 mjy beam -1 , while the contours are at 5, 10, 20, 30, 60, 90, 120 and 
150 mjy beam -1 . The FWHM of the Gaussian restoring beam is shown at the lower right of each image. The diagonal 
line above the SNR runs parallel to the Galactic Plane. Various regions of interest are marked. The ellipse in the left panel 
shows the lcr error ellipse for the pulsar's birth-place. 



of the RCW 89 region with the ROSA T Hig h Resolution 
Imager (HRI) data of Brazier & Becker (1997). The ring of 
knots seen at the radio peak of RCW 89 shows a marked 
correspondence with the X-ray emission. 

To the north of the SNR is source 1, an unresolved 
core surrounded by a patchy plateau of emission of diameter 
~3 arcmin. This source was classifie d as a SNR candidate, 
G320.6-00.9, by Whiteoak & Green ( |l996| ). Its position and 
20- and 6-cm flux densities are given in Table ^. Its flux den- 
sity at 6 cm is difficult to estimate, as the correction for the 
primary beam response is large. Furthermore, the extended 
component was not imaged by the 6-cm u — v coverage, and 
thus the corresponding flux density pertains only to the core 
region. 

4.2 Polarization 

An image of linear polarization L at 20 cm is shown in Fig- 
ure |(| showing diffuse patchy polarization over the SNR, 
with noticeable peaks in regions 4, 6, 8 and 9. Polarized 
intensity at 6 cm (not shown) shows similar peaks but at 
higher resolution. In Table ^| we summarise the fractional 
polarization in different regions, showing a general increase 
at 6 cm compared with 20 cm. 

ATCA continuum observations involve the recording of 
many adjacent frequency channels. These can be used to ex- 
tract the Faraday rotation across the observing band, from 
which an accurate rotation measure (RM) can be derived 
(Gaensler et al. 1998a). RMs determined from Faraday ro- 
tation across the 20-cm band are generally low (|-RM| < 



Figure 3. Grey-scale image of part of SNR G320.4-01.2 at 6 cm. 
The image has been tapered towards the edges to give uniform 
noise across the image, and shows emission in the range to 
20 mjy beam -1 . The FWHM of the Gaussian restoring beam is 
shown at lower right. 



100 rad m -2 ) except in regions 2 and 6 and on the south- 
ern edge of region 4. The RM of the SNR at the position of 
B1509-58 is 210±30 rad m -2 , very similar to the value of 
215 ± 2 rad m -2 for the pulsar (F. Crawford 1998, private 
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Figure 4. Radio/X-ray comparison of the region near PSR B1509— 58. The grey-scale corresponds to a sub-region of 
Figure^, but showing 20-cm emission in the range 1 to 10 mjy beam -1 . The contours represent the ROSAT PSPC data 
of Greiveldinger et al. flippy , with levels (in arbitrary units) at 1.5, 2, 3, 15, 30. Corrugations seen in the radio emission 
are low-level artifacts of the deconvolution process. 



Table 3. Peak linear polarization of regions indicated in Figure 



Region 


Fractional polarization (per cent) 




20 cm 


6 cm 


1 


< 3 


< 3 


2 


17 


40 


3 


5 


25 


4 


10 


18 


5 


5 


20 


6 


60 


120 


7 


13 


28 


8 


7 


30 


9 


32 





communication). The pulsar is too weak to contaminate the 
value measured for the remnant. 

The fractional bandwidth at 20 cm is insufficient to ac- 
curately determine intrinsic position angl es (s ee Gaensler 



et al. 1998a), while as discussed in Section 3.1, the 6-cm L 
image cannot be easily deconvolved. Thus we cannot derive 
intrinsic position angles of linear polarization from the SNR, 
and re fer th e reader to the lower resolution results of Milne 
et al. ( 1993[ ). In particular, we note that their results show 
that the magnetic field in the highly polarized region 6 is ori- 
ented along the axis defined by the collimated X-ray feature 
and the associated radio channel both seen in Figure [| 

The clumpy nature of polarized emission at 20 cm, along 
with the significant variations in RM seen across the rem- 
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Figure 5. A high-resolution radio/X-ray comparison of the 
RCW 89 region. The grey-scale delineates 6-cm ATCA obser- 
vations, while con tours correspond to the ROSAT HRI data of 
Brazier & Becker (1997), smoothed to the resolution of the 6-cm 



image. Contour levels (in arbitrary units) are at levels of 5, 10, 
15, . . ., 40. 
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Figure 6. Polarized emission from G320.4-01.2 at 20 cm. The 
grey-scale runs from 1.2 to 5.0 mjy beam . Contours represent 
total intensity at 20 cm smoothed to 50-arcmin resolution, and 
are at levels of 20, 100 and 200 mjy beam . 



nant, suggest that the low levels of L/I seen at 20 cm can 
be explained by beam depolarization, in which the position 
angle of polarized emission varies across the sky on scales 
smaller than the beam. As expected, the fractional polar- 
ization is significantly higher at 6 cm, a result of both the 
increased resolution and the reduced Faraday rotation at 
shorter wavelengths. Since our RMs are derived using up 
to 13 independent data points across the 20-cm band, it is 



unlikely that there is any ambiguity in their determination. 
There is reasonable agreement between the R Ms determined 
here and those derived by Milne et al. (1993) using only two 
data points, and thus one can have confidence in the intrin- 
sic position angles and rotation measures which Milne et al. 
derive. 

At 6 cm, radio emission near PSR B1509-58 (region 6) 
has the unphysical value of L/I — 1.2. This can be under- 
stood in terms of the spatial filtering applied by an inter- 
ferometer to the sky distribution of emission. In 6-cm total 
intensity, region 6 is poorly imaged by the ATCA because 
it is dominated by extended structure to which the tele- 
scope is not sensitive. But in polarization, variations in po- 
sition angle (either intrinsic to the source or resulting from 
differential Faraday rotation) shift power in the polarized 
images into smaller-scale structure to which the interfer- 
ometer is sensitive. Polarized emission thus appears to be 
of greater surface brightness than the total intensity (see 
Wieringa et al. 1993). Although the fractional polarization 
of region 6 is also high at 20 cm, it does not suffer from this 
effect as the observ atio ns are sensitive to all relevant spatial 
scales (see Section . The high fractional polarization at 
20 cm thus shows that there is minimal beam depolarization 
in this region, indicating a well-ordered magnetic field. 

No linear polarization is detected from G320.6-00.9 
(source 1 in Figure ^). In circular polarization, the 20-cm 
image shows emission only from PSR B1509-58, while the 
6-cm image contains no discernible sources. 

4.3 Spectral Index 

Appropriate spatial filtering of the 20- and 6-cm data results 
in a pair of images both of which contain emission on scales 
between 24 arcsec (the resolution at 20 cm) and 520 arcsec 
(corresponding to the shortest spacings at 6 cm), and have 
no structure outside this range. T-T plots were produced 
for the various regions of G320.4-01.2 marked in Figure [| 
Several plots were produced for each region by shifting the 
box in which data were considered by a few pixels each time. 
In all cases adjacent boxes produced consistent results. 

The results of the T-T fitting are shown in Table §. No 
value of a could be calculated for source 1 using this method, 
both because of the small number of pixels available for the 
fit and because of the large uncertainties in the primary 
beam correction near the edge of the field; simply using the 
total flux density at 20 and 6 cm we estimate a|° — —0.5 ± 
0.1. Region 6 was only barely imaged at 6 cm because of lack 
of short spacings in the array; we attempted to calculate a 
spectral index for it using our 20 cm d ata t ogether with 
the 36 cm image of Whiteoak & Green (1996), but its low 



surface brightness and the nearness of the two wavelengths 
prevents us from constraining its spectrum any more tightly 
than — 1 < «2o < 0- Region 9 was outside the region of the 
6-cm mosaic, and we roughly estimate a^o = —0.5 ± 0.2. 

We expect the flux density scale of the ATCA images 
to to be accurate to ~3 per cent. While uncertainties in ab- 
solute flux calibration will always limit the accuracy with 
which a itself can be calculated, this effect biases the spec- 
tral index of all regions by the same amount, and thus does 
not affect our ability to comment on spatial variations in a. 
Errors caused by applying a single primary beam correction 
on images formed from a large fractional bandwidth (as is 
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Table 4. Spectral indices for regions marked in Figure | 
Region 



20 
*6 



1 


(see text) 


2 


-0.40 ±0.03 


3 


-0.45 ±0.04 


1 


-0.52 ±0.04 


5 


-0.42 ±0.04 


6 


(see text) 


7 


-0.34 ±0.03 


8 


-0.33 ±0.03 


9 


(see text) 



the case at both 20 and 6 cm) can bias a towards the edges of 
the field, but the effect is negligible (Aa 0.01) compared 
to the uncertainties quoted. 

Thus results from Table ^ indicate that the spectrum 
of region 4 is 2-3a steeper than that of regions 2 and 5, and 
3-4cr steeper than the spectra of regions 7 and 8. 



4.4 Hi line 

One usually compares absorption spectra with emission 
spectra from an adjacent patch of the sky. However, a use- 
ful emission spectrum is difficult to extract from our ATCA 
data because of insufficient u — v coverage on the large spa- 
tial scales where emission dominates. We thus compare our 
results to the emis sion s pectrum in this direction obtained 
by Caswell et al. ( 1975 ). This spectrum shows significant 
emission at LSR velocities of -70 km s - , -55 km s , from 
-15 to ±30 km s -1 , then at ±80 km s -1 . 

Of unrelated sources near G320.4-01.2, only G320.6- 
00.9 (source 1 in Figure []) has enough signal-to-noise ra- 
tio to obtain a useful absorption spectrum. This spectrum 
is shown in Figure |t], where significant absorption is seen 
at -75 and -55 km s -1 and between -10 and ±20 km s , 
with a lower level feature at ±65 km s _1 . We thus put 
a lower limit on the systemic velocity of G320.6-00.9 of 
V L = ±20 km s -1 . We use the best fitting model for Galac- 
tic rotation of Fi ch, Blitz fc Stark (1989), a dopt standard 
IAU parameters (Kerr & Lynden-Bell 1986) for the solar 
orbital velocity (Oo = 220 km s - ) and distance to the 
Galactic Centre (Ro = 8.5 kpc), and assume an uncertainty 



of ±7 km s in systemic ve locities, representa t ive of the 



random motion of H I clouds ( B haver et al. 1982j ; Belfort & 
Crovisier 1984). We can thus derive a lower limit on the 
distance to G320. 6-00.9 of 14.8 ± 0.7 kpc. 

Thus G320.6-00.9 is at a large distance, is unpolarized, 
and has a comparatively steep spectral index, and its core is 
unresolved at our highest resolution (2 arcsec using the 6-km 
antenna at 6 cm). It therefore seems likely that this sourc e 
is not a SNR as proposed by Whiteoak & Green (1996), 
but rather is a background radio galaxy. The nature of its 
extended envelope is unclear. 

In Figure ^ we also show H I spectra towards three parts 
of G320.4-01.2. Regions 2 and 4, both in the northern part 
of the SNR, show the same main features in absorption: 
strong absorption at -55 km s -1 and then at -15 km s _1 
(region 2) or -5 km s -1 (region 4). Other spectra towards the 
northern part of the SNR are similar, all showing absorption 
out to -55 km s . Region 8, in the south of the SNR, shows 



significant absorption only at -25 km s -1 and -5 km s . 
A feature at -60 km s~ is at the limit of significance. The 
signal-to-noise of absorption towards other parts of the SNR 
was too low to make any useful measurements. 



4.5 OH and H140a 

Neither OH nor H140a emission was detected towards 
G320.4-01.2. The 3-<r upper limit on the surface brightness 
of OH is 23 mjy beam -1 at 4 km s -1 resolution, consistent 



with the survey results of Green et al. (1997) but at a slightly 
more stringent limit. Limits for H140q are 6 mjy beam - 
(15 km s -1 resolution) and 4.5 mjy beam -1 (40 km s -1 res- 
olution). At the peak of the radio emission from RCW 89, 
the H140a non-detection corresponds to an upper limit on 
the line-to-continuum ratio of ~4 per cent. 



5 DISCUSSION 

5.1 How many SNRs? 

Although there are several distinct components to the radio 
emission from G320.4-01.2, each part can be matched up to 
part of the X-ray emission in the region. 

(i) Based on the correspondence shown in Figure ^, the 
north-western radio region (which appears to be a single 
object based on the similarity of Hi spectra for regions 2 
and 4) can be associated with the X-ray region "NN" . 

(ii) The morphological correspondence seen in Figure ^ 
between radio region 6 and the X-ray region "SN" argues a 
connection between them. 

(iii) Figure ^ demonstrates a match between the "T" re- 
gion in X-rays and regions 7 and 8 in the radio. 



Trussoni et al. (1996) have shown that the spectra of 
and absortion towards these X-ray regions are consistent 
with their all being part of a single SNR. Since each of the 
main radio regions corresponds to part of this X-ray rem- 
nant, we therefore argue that the entire G320. 4~01.2 region 
in radio and X-rays is a single SNR. 



5.2 A distance to G320.4-01.2 

Absorption spectra for regions 2 a nd 4 of the SNR are con- 
sistent with that of Caswell et al. (197E): absorption is seen 
at -55 km s -1 but not at more negative velocities. We thus 
adopt Vl = —55 km s , corresponding to a lower limit on 
the distance of 3.8 ± 0.5 kp c. An upper limit i s the tan- 
gent point at ~ —70 km s -1 (Caswell et al. 1975 



Kerr et al 



1986), which is at a kinematic distance 6.6±1.4 kpc. Region 
8's spectrum can only be considered consistent with those 
of regions 2 and 4 if absorption at -60 km s -1 is genuine. 
However, region 8 is several arcmin away from the other 
sources, and the distribution of H I clouds along the line of 
sight may be quite different. Indeed, values for the hydro- 
gen column density towards the SNR inferred from photo- 
electri c absorption indicat e significant gradients across the 
SNR ( Trussoni et al. 1996 ). Thus the Hi spectra can all be 
considered consistent with a distance for SNR G320.4-01.2 
of 5.2 ± 1.4 kp c. Th is agrees with the distance derived by 
Caswell et al. (1975), but has more realistic uncertainties. 
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Figure 7. H I absorption spectra towards four regions in the field marked in Figure | 
at the 6<r level, where a is the noise in line-free channels. 



The dashed line marks absorption 



In further discussion we adopt a distance to the SNR of 
5ds kpc. 



thors to propose that the pulsar's wind is focused into rel- 
ativistic jets or outflows, generated al ong an axis aligned 
approximately north-west /south-east (Seward ct al. 1983 



5.3 



radio emission near Jr"SJrt hsitjuy t>» 



Manchester fc Durdin 1983; 



Region 6 has the same RM as the pulsar (quite different 
from that of RCW 89), has a distinctly higher fractional 
polarization than RCW 89, and shows a morphological cor- 
respondence with X-ray emission associated with the pulsar 
nebula. Thus there is good evidence that region 6 is directly 
associated with PSR B 1509- 58, as originally suggested by 
Manchester & Durdin (198S). Despite this claim, we would 
argue that this region of radio emission is not a radio PWN 
in the traditional sense, since the radio emission is not cen- 
tred on the pulsar, nor are the X-ray-bright regions bright 
in the radio. Indeed any radio PWN is expected to be faint: 
comparison with the Crab Nebu la gives an expected flux 
density at 20 cm of ~100 mjy (Brinivasan, Bhattacharya 
&i Dwarakanath 1984), while extrapolation from the X 



s pectrum nf ttic segjau i™pn»a a fl " v density ~ j1 ™ T y 

(gewarl et al. 19S4[h If we a.ssume that, this radio PWN is of 



the same extent as the X-ray PWN (although in fact we ex- 
pect the radio source to be larger because of the longer syn- 
chrotron lifetimes), we expect a 20-cm mean surface bright- 
ness at 1-arcmin resolution of 0.02-2 mjy beam - , com- 
pletely undetectable against the background in this region 
of ^>20 mjy beam -1 . 

The morphology of the X-ray PWN, together with the 
appearance of the surrounding SNR, have led several au- 



Tamura et al. 199C 



Brazier & 



Becker 1997). There are collimated X-ray features along this 
axis on both sides of the pulsar. In the south-east, Figure ^ 
shows that radio emission is enhanced along the sides of the 
collimated X-ray feature, but not within it. We argue that 
this anti-correspondence can be explained if the X-ray fea- 
ture is indeed a jet or outflow along the previously proposed 
axis; radio emission is then produced by shocks generated 
in a cylindrical sheath around the jet, as seen for SS 433 
(Hjellming & John ston 1986, 1988) and possibly for the X- 
ray jet in Vela X (Frail et al. 1997). A jet interpretation is 



also favoured by the polarization properties of the sheath: 
the magnetic field in the region is well-ordered and orien ted 
parallel to the axis of the proposed outflow (Section 4.2 
above; Milne et al. 1993), similar to the field structu re seen 
surroundin g the jets associated with both the Crab (Hester 



et al. 1995) and Vela (Markwardt fc Ogelman 1995 



Milne 



1995 ) pulsars. 

We assume the radio sheath to be a hollow cylinder, 
the inner and outer diameters of which correspond to the 
diameters of the X-ray jet and radio sheath respectively. 
From ROSAT PSPC and HRI data we estimate the inner 
diameter of the sheath to be 60-75 arcsec, while from 20-cm 
ATCA data we determine an outer diameter of 210 arcsec. 
We assume that the jet is inclined to the line of sight at 70° 
( Brazier fc Becker 1997 ) and that the radio emission in the 
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sheath is optically thin. We then expect the radio emission 
along lines of sight through the X-ray jet to be 0.7 as bright 
as the radio emission along the edges of the jet, consistent 
with the observed value of 0.5-0.8. The sheath is too faint for 
us to put any useful constraint on its spectral index, but we 
expect it to have a —0.5, indicative of shock acceleration 
(iFrail et al. 19971). 



The abrupt fading of the sheath to the south-east is 
difficult to explain, but appears to coincide with a similar 
termination in X-rays seen in Figure ^. Clearly, conditions 
within or surrounding the jet change in some way beyond 
this point. 

The cross-shaped morphology of the X-ray PWN and 



its rese mblance to t he Crab Nebula ( Hester et al. 1995 ; Bra- 
zicr & ftecker 1997| ) both argue for the presence of a second 
jet to the north-west of the pulsar. This jet appears to be 
less collimated than the south-eastern jet and has no obvious 
radio depression associated with it. While there is no clear 
explanation for this difference, we note that conditions on 
the north and south sides of the pulsar may be quite differ- 
ent (see Section 5.7 below), and that asymmetries are also 



observed in the outflows associated with the Crab and Vela 



pulsars (Pelling et al. 1987; Markwardt fc Ogelman 1995) 



5.4 The knots in RCW 89 



The north-western jet which we have just discussed appears 
t o map onto the ring of X-ray /radio knots seen in Figure ^ 
(Brazier & Becker 1997). Thus rather than argue that the 
entire RCW 89 region is p o wered by the pulsa r outflow 



(Manchester & Durdin 1983; Tamura et al. 1996), we now 



consider the alternative that the interaction region is con- 
fined just to the ring of knots. 

The radio knots within RCW 89 are significantly lin- 
early polarized, show no recombination lines and have the 
steepest spectral index of any emission from the SNR. Thus 
there seems little doubt that the emission mechanism associ- 
ated with the radio features is synchrotron radiation. Taking 
the brightest knot to be a sphere of radius 15 arcsec and flux 
density 0.2 Jy at 20 cm, we can u se standard minimum en- 
ergy arguments ( Pacholczyk 197C| ) to infer a magnetic field 

/iG, where is the filling 



2 / 7 rh~ 2 / 7 dt 11 '' 



for it of 1.5(1 + ky> 7 <f> 
factor of emitting particles and fields and k is the ratio of en- 
ergy in heavy particles to that in electrons. This is not signif- 
icantly different from the field stren gths of 7 /iC and 2. 5 /iG 
inferred for the ov erall X-ray PWN ( seward et al. 1984) and 
for the X-ray jet (Tamura et al. 1996) respectively. The cor- 



responding total energy of this knot in magnetic fields and 
relativistic particles is 1.3(1 + k) 4/7 4> 3/7 d 17/7 x 10 42 erg, a 
tiny fraction of the total spin-down energy released by the 
pulsar in its life-time (> 10 48 erg). The steep spectrum of 
the knots is the opposite to what we might expect if injec- 
tion of particles from the pulsar was directly contributing 
to their emission (cf. Frail et al. 1994b). Rather, it seems 
that the collision between the jet and pre-existing material 
drives shocks into the latter at the working surface, produc- 
ing clumps of enhanced synchrotron emission. 

X-ra y observations of RCW 89 as a whole sho w it to be 
thermal ( Reward et al. 1984 ; Tamura et al. 1996), but the 
X-ray/radio correspondence seen in Figure ^ suggests that 
X-ray emission from the knots might be synchrotron emis- 



sion as in the radio. A single power law of spectral index 
—0.6 Ss a 5; —0.5 (from Table ^) can indeed be extrapolated 
between the radio flux of the knots as observed he re an d the 
X-ray count rates quoted by Brazier & Becker (1997). No 
spectral break as might be caused by synchtrotron losses is 
required, which implies a magnetic field strength < 10 /iG 
consistent with the determination from minimum energy 
above. While RCW 89 is predominantly thermal in X-rays, 
the available spectral data lack the spatial resolution to de- 
termine whether the thermal spectrum of the region corre- 
sponds to the knots as well as the more extended component 
in which they are embedded. C omparison of the coun t-rate 
for the entire RCW 89 region (Trussoni ct al. 199Q with 
those of the knots alone ( Brazier fcFjecker 1997|) shows that 
soft X-ray emission is dominated by the extended compo- 
nent, and thus at low energies a power-law spectrum associ- 
ated with th e knots might be hid den. Indeed observations in 
hard X-rays (Tamura et al. 1996) show non-thermal emission 
at their position. However, the spatial resolution of these 
data is not sufficient to separate the knots from more dif- 
fuse emission in the area; high-resolution observations with 
AXAF will be required to clarify the issue. Whether thermal 
or non-thermal, the energy in the knots is far smaller than 
that in the whole RCW 89 region, and so can be comfortably 
accounted for by the pulsar spin-down. 

The horseshoe-shaped plateau in RCW 89 appears to 
have a similar shape to the ring it envelopes, and may rep- 
resent material diffusing away from the point(s) of impact. 
However, the size of this region is then somewhat large: from 
the ring to the edge of this plateau region is 2.5 arcmin, 
which requires emission to have travelled outwards from the 
ring at a projected mean velocity of 2100ds km s _1 . One 
possibility suggested by Tamura et al. (1996) is that the 



pulsar jet has precessed or narrowed in opening angle. 

To summarise, the problems raised at the end of Sec- 
tion 1.2 can be resolved if just the X-ray/radio knots, rather 
than the whole of RCW 89, are interpreted as the point of 
interaction between the pulsar outflow and the SNR. The en- 
ergy in these knots can be easily accomodated by the pulsar 
spin-down, and it is easier to understand how these features, 
rather than the whole extended nebula, might be the termi- 
nation of an energetic outflow. A more likely possibility to 
account for the surrounding RCW 89 region is that it orig- 
inates from the SNR blast wave, and is not related to any 
outflow from the pulsar. 



5.5 The south-eastern component of the SNR 

It can be seen in Figure |^ that regions 8 and 9 of the SNR 
contain compact polarized clumps, resembling the knots 
seen in RCW 89. These features may be evidence for an 
interaction between the south-eastern jet and the south- 
eastern half of the SNR, although, compared to that with the 
RCW 89 region, the impact is greatly reduced. In Figure ^ 
we show an e xtension of the model proposed by Brazier & 
Becker (1997), where we speculatively project an additional 



cone of emission onto the south-eastern part of the SNR. 



t This] 

Durdin (1983) and Manchester (1987), but maps onto a smaller 



oiect ion is qualitatively simil ar to that of Manchester & 



region of each half of the SNR. 
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Figure 8. Speculative geometry for outflow in G320.4— 01.2. Con- 
tours represent 20-cm continuumm, and are at 5, 20, 60, 120 and 
150 mjy beam . 



It is interesting to note that this second cone happens to 
intersect the SNR at the location of these two highly polar- 
ized clumps. It can be seen in Figure |l]that while the X-ray 
emission joining the "SN" and "T" regions might make up 
one edge of this cone, the other edge is missing. However, 
the latter region lies right along a rib of the PSPC and so 
may have had minimal exposure. 



5.6 Spectral index variations 



Results presented in Section 4.3 indicate that the compact 



knots at the peak of RCW 89 have a steeper radio spec- 
trum than the diffuse emission surrounding them, which 
are again steeper than the south-eastern part of the SNR. 
While the phenomenon of spatial spectral index variations 
in SNRs is not well understood, there is the suggestion that 
young remnants generally show a stee per spectrum in their 



brigh test, most compact components (Anderson & Rudnick 



1993), as observed here. For example, in G260.4-03.4 (Pup- 



pis A) the steepest spectrum part of the shell is where it has 
been apparently dec elerated by a collision with a dense cloud 
(bubner et al. 1991), while in Gil 1.7-02.1 (Cassiopeia A), 



the steep spectrum regions tend to be bow-shocks drive n 
by clumps of fast-moving material (Anderson et al. 1991). 



We note that both these effects may be occurring in G320.4- 
01.2: si gnificant deceleration of the RCW 89 region is argued 



for in S pction 5.7| below, while the pulsar outflow may accel- 
erate small clumps of material to high velocity, which then 
drive bow-shocks as in Cas A. 



5.7 Finale: what is SNR G320.4-01.2? 

The dis tance to the SNR of 4.2 kpc as derived by Caswell 
et al. ( 1975j ) was used in the calibration of the electron 
distribution model of Taylor & Cordes (1993). Thus com- 
parisons between the pulsar's distance as derived from this 
model and any distance estimate to the SNR are not en- 
tirely meaningful. Rather, we estimate the pulsar's distance 
as follows. 

While PSR B1509-58 itself is too weak to show Hi ab- 
sorption, PSRs B1240-64, B1323-58, B1323-62, B1356-60 
and B1557-50 are all within 20° of PSR B1509-58, have 
comparable dispersion measures and have independent dis- 
tance estimates from H I absorption ( Frail & Weisberg 199C ; 

1996D- Distances and dispersion measures 



Saravanan et al. 



for these pulsars suggest a mean electron density along this 
line of sight of between 0.03 and 0.06 cm -3 , corresponding 
to a distance for PSR B1509-58 in the range 4.2 to 8.4 kpc. 

This distance for the pulsar is consistent with the dis- 
tance we have determined for SNR G320.4-01.2 through 
Hi absorption, resolving any previous discrepancy. Further- 
more, despite the large variation in RM across the SNR (- 
200 to +400 rad m~ 2 ), its RM at the position of the pulsar 
agrees with that of the pulsar itself. Finally, we discussed in 



Section 5.4 that the pulsar appears to be interacting with the 
surrounding SNR. Thus we argue strongly that PSR B1509- 
58 and SNR G320. 4-01.2 are associated. The age discrep- 
ancy outlined in Section [l^ must then be resolved. 

A first possibility to account for this apparent age differ- 
ence is that B 1509-58 has had a n unusual spin-down history , 
and is much older than it seems (Blandford & Romani 1988). 



However, all observations suggest that B 1509-58 is a typi- 
cal young pulsar, but that G320.4-01.2 is an anything-but- 
typical remnant. Although the possibility that B1509-58 is 
old cannot be ruled out, disagreement between a relatively 
well-determined pulsar age and a crudely estimated SNR age 
gives little cause to look to the pulsar as the culprit. Thus 
in future discussion we assume that the age of the system is 
5sl700 yr. We now attempt to account for this age in terms 
of the properties of the supernova and its environment. 

We take the south-eastern component of the SNR (in- 
corporating regions 7, 8 and 9) as that part most represen- 
tative of a "normal" shell SNR. Assuming the site of the 
supernova explosion to be within the 1-cr error-box derived 
from the pulsar's proper motion, we estimate the radius of 
the south-eastern component of the SNR to be 17 ± 2 ar- 
cmin. For an SNR age of 1700 yr, this implies a mean ex- 
pansion velocity of (14 ± 2)d 5 x 10 3 km s" 1 . If this com- 
ponent represents freely expanding material, we find a ki- 
netic energy for the explosion of E51 — (2.0 ± 0.3)M c jd§, 
where M e j is the ejected mass in units of Mq and -E51 is 
in units of 10 51 erg. Alternatively, if we assume expansion 
in the adiabatic (Sedov- Taylor) phase then the size of the 
SNR implies (Esi/m) ~ (1100 ± 600)4, where m is the 
ambient density into which this component of the SNR is 
expanding. To fully enter t his phase the SNR must sweep 
up 20 times its own mass ( Fabian, Brinkmann fc Stewart 



198S ; Dohm- Palmer & Jones 1996), implying a density ni > 
(0.012 !± 0.003) Afcjd" 3 cm" 3 to give E 51 £ (3.3 ± 2.9)M ej d|. 

Thus regardless of the evolutionary state of the SNR, 
an age of 1700 yr implies a value of Esi/Md which is 
10-100 times higher than typical values ( Berkhuijscn 198*^ ; 
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Smith 1988). This problem has been overcome in the past 
by invoking a large [Es\ ~ 10 — 100) kinet ic energy for 



- Esi/M cj ~ 2 as argued for above. Scaling to Figure 3 of Gull 
( 1973| ), we find that the north-western component of the 



the supernova ( |Seward et al. 1983| |Katz 1983[ ). While such SNR is well into the adiabatic phase, and that the ambient 

density in this region is in the range ni — 1 — 5 cm -3 . This 



a value for E$i seems uncomfortably high, it is compara- 
ble with (model-dependent) estimates for the e nergy re- 
leased i n the recent Type Ic supernova SN 1998bw (Iwamoto 
et al. 1398; Woosley, Eastman & Schmidt 1998). However, 



we propose the more likely explanation that the progeni- 
tor of PSR B1509-58 was a helium star. Such a star ex- 
plodes in a Type lb supernova, ejecting a low amount of 
mass (M e j ~ 1 .5) but with a typical supernova energy 
(Ebi ~ 1 — 2) (Woosley, Langer & Weaver 1995), as re- 



quired here. While the optical fila ments in RCW 89 show 
no evidenc e of helium rich eje cta ( Dopita, Mathewson & 



is consistent with the non-detection of OH, as the 1720 MHz 
maser line will b e shock-excit ed into emission only at much 
higher densities ( Elitzur 197£| ). 

To summarise, we find that a number of observed prop- 
erties of the SNR, particularly its bilateral appearance, faint 
radio PWN, asymmetric distribution about the pulsar and 
large size (and possibly also the slow optical expansion of 
filaments in RCW 89) can all be explained in a model in 
which an explosion of high energy or low mass occurred 
near the edge of a cavity elongated parallel to the Galac- 
Ford 1977; Seward et al. 1983), we show below that this tic Plane. The density of the cavity is ~ 0.01 cm -3 : the 

south-eastern component of the SNR has expanded rapidly 
across the cavity and has recently collided with the other 
side, while to the north-west, the shock has encountered 
denser (1-5 cm -3 ) material, causing the SNR to decelerate 
significantly. To the north-east and south-west the shock is 
still propagating through the cavity and produces no observ- 
able emission. 

While at low resolution G320.6-01.6 appears to be an 
extension of G320.4-01.2 (Milne et al. 1993), it is impossi- 
bly large to have expanded to its current size in 1700 yr. 
Its morphology and surface brightness are in complete con- 
trast to those of G320.4-01.2, and we conclude that it is an 
unrelated older SNR along the same line of sight. 



component of the SNR is considerably evolved and is now 
dominated by swept-up material. Spectroscopy of any opti- 
cal emission which can be associated with the south-eastern 
component of G320.4-01.2 would be of great interest. 

Helium stars start off as extremely massive stars, but 
lose their hydrogen envelopes either through a strong stel- 
lar wind or through transer of mass to a binary companion 
(Woosley, Langer & Weaver 1993, 1995). The latter possi- 
bility is interesting, particularly a s the O star Muzzio 10 
(Muzzio 1979; Seward et al. 1982 ) is just 18 arcsec north 
of the pulsar, and is at a comparable distance (~5 kpc; 
Arendt 1991). We thus speculate that Muzzio 10 is the for- 
mer companion of PSR B1509-58. This possibility is easily 
tested: assuming a pulsar age of 1700 yr, the association pre- 
dicts a proper motion for the pulsar of 11 mas yr -1 along 
a position angle 168° (N through E), corresponding to a 
transverse velocity 260ds km s~ . 

It is clear from the radio data that there is no radio 
emission which can trace an entire shell and thus join the 
two components of the SNR. Gaensler (1995) has argued 
that the "bilateral" or "barrel" morphology which results 
(and the observed alignment of the bilateral axis with the 
Galactic Plane) is produced when a SNR expands into an 
elongated low-density cavity, radio emission being produced 
where the shock has recently encountered the walls. Thus a 
bilateral appearance can be considered evidence that ni is 
low (~ 0.01 cm -3 ); such a value f or ni can also accou nt for 
the faintness of the radio PWN ( Bhattacharya 199C), and 
is consistent with the cavities expected due to the presence 
of other massive stars in the region (Lortet et al. 1987). For 
Aicj = 1.5, this ambient density corresponds to the south- 
eastern component beginning to enter the adiabatic phase. 

Limits on the pulsar proper motion as indicated in Fig- 
ure ^ demonstrate that the SNR is distributed highly asym- 
metrically about the site of the explosion. Intuitively, this 
suggests that the supernova explosion was near the edge of 
the proposed cavity, and that, while the south-eastern half 
of the SNR has expanded relatively unimpeded, the north- 
western half has decelerated significantly. A high-density en- 
vironment for this component is also argued for by its high 
X-ray and radio brightness and by the small prope r motion 



of associated optical filaments (Seward et al. 1983 
ester 19^3). 



Manch 



Gull (1972) has modelled the deceleration of a young 
SNR, and provides means of scaling his solution to arbitrary 
values of density and Esi/M e $. We take the radius of the 
north-western component to be 7 ± 2 arcmin, and adopt 



6 CONCLUSION 

We have reported on an extensive set of ATCA observations 
of G320.4-01.2 and PSR B1509-58. The main results of this 
study are as follows. 

(i) The disparate radio components of G320.4-01.2 are 
all part of a single SNR at a distance of 5.2 ± 1.4 kpc and 
with an age of ~1700 yr. 

(ii) PSR B1509-58 is physically associated with G320.4- 
01.2. The pulsar emits twin jets or collimated outflows of 
relativistic particles, one of which interacts with the SNR in 
the form of radio/X-ray knots within RCW 89. 

(iii) SNR G320.4-01.2 was formed in a supernova of high 
kinetic energy or low ejected mass (Esi/M e j ~ 2) which 
occurred near the edge of a low-density cavity. 

A variety of further observations can support or refute 
these conclusions. For example, the cavity we propose may 
be visible in H I emission, while forthcoming AX AF obser- 
vations will allow a detailed study of the physical conditions 
within the knots. 

While we have offered answers to the questions out- 
lined in Section [lj the picture is far from complete. First, 
we lack an understanding of the physical details of how the 
pulsar outflow is generating the peculiar knots in RCW 89, 
and of how the pulsar is interacting with the south-eastern 
half of the SNR. Similar interactions through outfl ows or 
lets have been proposed for a variety of other SNRs ( Roger 
et al. 19851 ; | Vlanchester 1987] ; |Gaensler, Green fc Manchester] 
1998b), although an associated pulsar is yet to be detected 
in most cases. We have also sidestepped the whole issue 
of why PSR B 1509-58 generates such large-scale jets and 
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how general this situation might be. Certainly the theoreti 
cal expectation has been for s ome time that p ulsars should 
generate collimated outflow s (Benford 1984 



Michel 1985 



Sulkanen fc Lovelace 199C| ). The Crab, Vela, PSR B1509- 
58 and possibly PSR B1951+32 (Hester 1998; J. J. Hester 
et al. in preparation) all show evidence that this is indeed 
the case, while various jet-like fea tures have been proposed 
around other pulsars (iBell 1997). We must now therefore 



accept that spherically symmetric pulsar winds are a gross 
over-simplification . 
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